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Abstract 



We study the dynamics of a Brownian quantum particle hopping on an infinite lattice with a spin degree 
of freedom. This particle is coupled to free boson gases via a translation- invariant Hamiltonian which is linear 
in the creation and annihilation operators of the bosons. We derive the time evolution of the reduced density 
matrix of the particle in the van Hove limit in which we also rescale the hopping rate. This corresponds to 
a situation in which both the system-bath interactions and the hopping between neighboring sites are small 
" and they are effective on the same time scale. The reduced evolution is given by a translation-invariant 

, Lindblad master equation which is derived explicitly. 

KEY WORDS: out-of-equilibrium quantum statistical physics, open quantum systems, weak coupling limit, 
singular coupling limit, quantum Brownian motion. 



1 Introduction 



- , — i ■ 

The irreversible dynamics of a quantum system coupled to infinite baths is often described by determining 
the time evolution of the reduced density matrix of the system, the latter being obtained by tracing out the 
bath degrees of freedom in the system + bath state. Under certain approximations (including a Born-Markov 
approximation), this density matrix is the solution of a Lindblad master equation [TJH]. We are aware of three 
mathematically well-defined ways to derive such a Lindblad equation starting from the Hamiltonian dynamics 
of the system and baths [3J HI [5]: the weak coupling limit, the singular coupling limit, and the low density 
limit. The weak coupling limit goes back to [5] and it was put on a rigorous footing in a series of papers by 
E.B. Davies [7J[S]. It consists in letting the system-bath coupling constant A going to zero and rescaling time 
like t = A~ 2 t, with t > fixed. This limit enforces the separation of time scales 

t s « t R , t B < t R , (1) 

where ts (sometimes called the "Heisenberg time" ) is the time scale on which the system evolves in the absence 
of coupling with the baths, ts is the correlation time of the baths, and t R ss t is the time scale at which we 
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describe the dynamics, that is, the time scale on which the system evolves under the coupling with the baths 
(the "relaxation time" for systems converging to stationary states). The first time scale separation in ([TJ allows 
to perform the rotating wave approximation. The second time scale separation allows for the Born-Markov 
approximation. The singular coupling limit is the limit of delta-correlated baths, corresponding to 

«B « *S , «B « • (2) 

Such a limit, which is a quantum analog of the white noise limit for classical stochastic processes, has been 
analyzed rigorously in [9] 1 1 01 lllj . It is physically meaningful in the limit of large bath temperature. Finally, we 
refer the reader to [T^] for a description of the low-density limit. 

The weak coupling, singular coupling, and low-density limits have been applied and defined primarily for 
confined systems (typically atoms) coupled to free fermion or free boson baths. There is a compelling reason 
for this in the case of the weak coupling limit: the Hamiltonian of a confined system has discrete spectrum, 
and therefore a well-defined time scale ts (given by the minimum of the inverse level spacings); in contrast, 
extended systems may have continuous spectra, corresponding to arbitrarily slow processes in the uncoupled 
system dynamics, thus invalidating (TTJ (see, however, [T3] for a different approach). A physical example of 
this is diffusion, where the relevant time scale is set by a spatial scale. In contrast, the singular coupling limit 
remains well-defined for extended systems, as one can guess from inspection of and as we will illustrate 
in this article. Note that in the physics literature the dynamics of systems with arbitrarily large ts are often 
described by a Bloch-Rcdficld master equation with a time-dependent generator which is not of the Lindblad 
form (see [14] and references therein). This equation is perturbative in the system-bath coupling but does not 
include a rotating wave approximation. An explicit bound independent of ts on the difference between its 
solution and the exact reduced density matrix has has been obtained recently [15] . but otherwise this equation 
has, to our knowledge, not been analyzed rigorously. 

The derivation of the reduced dynamics of extended quantum systems is considerably more involved. In 
[TS], an extended system is studied in the scaling limit t = A~ 2 t, x = A~ 2 £, A — > (with r > and (el 11 
fixed) in which both the time t and the position x are rescaled. The scaling of space is dictated by the scaling 
of time, since on the microscopic scale the particle moves a distance of order A~ 2 in a time of order A~ 2 . In this 
limit, the resulting equation is a linear Boltzmann equation for the Wigner distribution of the particle. This 
framework has been extended to describe decoherence in position space in [17] and an essentially analogous 
result, with the weak coupling limit replaced by the low-density limit, was obtained in [18j . 

Let us also note that quantum systems coupled to infinite baths have been studied in the past fifteen year 
from another perspective. This approach, due to Jaksic and Pillet [19], uses operator algebras and spectral 
analysis to describe the dynamics of the system and bath at large time and small but finite coupling constant 
A. For confined systems, we refer the reader to the lecture notes collected in [20] and the references therein. 
Extended systems have been analyzed recently from a similar perspective in |21[ 122] . 

In this work, we consider an extended system coupled to bosonic baths. We are interested in the dynamics of 
the reduced density matrix of this system at long times and for weak couplings. The system is a quantum particle 
moving on an infinite lattice Z d , which has some internal degrees of freedom acting on a finite-dimensional Hilbert 
space. In the simplest case, the Hamiltonian of the particle is the sum of the discrete Laplacian — A on £ 2 (Z d ) and 
of a self-adjoint Hamiltonian S describing the internal degrees of freedom. One may also think of more general 
Hamiltonians coupling the position and internal degrees of freedom. The particle interacts with free boson 
gases via a translation-invariant Hamiltonian, assumed to be linear in the creation and annihilation operators 
of the bosons. We consider the following scaling limit: (i) the time is rescaled as t = A~ 2 r, where A is the 
particle-boson coupling constant and r > is fixed (ii) the particle Hamiltonian is rescaled as Hp = — A 2 A + S 
(iii) one takes the limit A — > 0. This scaling combines the weak coupling and the singular coupling limits: if the 
translational degrees of freedom are frozen, it reduces to the weak coupling limit for the internal state, whereas if 
one ignores the internal degrees of freedom it amounts to a singular coupling limit for the motional state, as will 
be explained below (see [3] [11]). In the latter case, however, the master equation is trivial. Indeed, by energy 
conservation the particle can only absorb or emit bosons with a vanishing frequency in the limit A — > 0; since 
such bosons have also a vanishing momentum and the total momentum is conserved, one has no momentum 
transfers between the particle and the baths; thus the coupling to the baths has no effect on the particle, except 
possibly for decoherence in the momentum basis. This is in fact the main reason why we consider a particle 
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with internal degrees of freedom even though we are primarily concerned with its motion on the lattice. Note 
that since the hopping strength is of order A 2 , we do not need a space rescaling as in [TBI H71 U5] . 

The above model allows for a tractable rigorous analysis in spite of the fact that we deal with a spatially 
extended system. Our main result states that, in dimension d > 2, the reduced density matrix of the particle 
converges in the aforementioned scaling limit to the solution of a Lindblad master equation which is determined 
explicitly. This equation contains the physics of dissipative extended systems, in particular diffusion (whose 
analysis is, however, not treated in this work, we refer the reader to [3TJ 123 for results in this direction). Its 
derivation requires much less mathematical complications than in the works [TBI [T71 [T8] . 

Some related models have been studied in [23j HH H3 US] . In particular, Vacchini and coworkers considered 
in a series of non rigorous works |24i I25| similar models but in a low density limit; they argue that the evolution 
of the particle density matrix (that is, not only of the associated Wigner transform) is governed in this limit by 
a Lindblad master equation. In [26j . drift and diffusion of an electron moving on a one-dimensional lattice and 
submitted to a static electric field has been studied in a model in which the coupling to the bath is simulated by 
repeated interactions with two level systems. Finally, we point out that our model can be viewed as a continuous 
version of a dissipative quantum walk [27]. 

The paper is organized as follows. We introduce the model in Section^ first at finite volume for the particle 
and baths and then by considering the infinite volume limit. Our results are presented and discussed in Section 
[21 together with two important examples. The last section 0] contains the proofs and some technical results are 
proven in the appendix. 

2 The model 

2.1 The quantum particle 

Our model consists of a quantum particle on the lattice 7L d coupled to free boson fields. In this subsection 
and the three following ones, we describe the model at finite volume^. We thus restrict the lattice to a finite 
hypercube with periodic boundary conditions and consider A = Z d /(2LZ) d with 1 < L < oo, d being the space 
dimension. We will often identify A with ] — L, L] d n Z d C 1i d . The infinite volume limit L — » oo will be taken 
in Subsection 12.51 The particle has translational degrees of freedom x £ A and an internal degree of freedom 
s = 1, . . . , N < oo, which may correspond to a spin or to an internal state of an atom or a molecule. The Hilbert 
space of the particle, Hp = £ 2 (A) ® C*, has finite dimension |A| x N, where |A| = (2L) d is the cardinality of 
A. The particle Hamiltonian Hp consists of a hopping term acting on £ 2 (A) plus a term governing the internal 
dynamics given by a self-adjoint operator S acting on C^, 

Hp = X a H^ op + S (3) 

where we identified S with (g) S. We have introduced in front of the hopping term a small parameter A Q 

playing the role of a hopping strength or of an inverse mass; A > will be chosen below to be the particle-boson 
coupling constant and a is a positive scaling exponent. Our most interesting result will correspond to a = 2 
and d > 2. Note that for a = oo and A < 1, i.e., A" = 0, the translation degrees of freedom can be dropped 
altogether and the dynamics takes place on C N . In the simplest setup, H^ op is (up to a minus sign) the discrete 
Laplacian on £ 2 (A), 

#hop = -A = - (\x){y\-\x){x\)®l c » (4) 

x,yeA,\x— s|a=1 

where {\x); x £ A} is the canonical basis of £ 2 (A), the Dirac notation \x)(y\ refers to the operator -0 >-> ">P(y)\x) 
from £ 2 {A) to £ 2 (A), and \x — y\\ — J2i=i mm fcez \x% — y% + 2kL\, i.e., we use periodic boundary conditions. 
In a more elaborate setup, H^ op will be modified such that the propagation of the particle may couple to the 
internal state. An Hamiltonian that accommodates this idea is presented in Section [3.31 

1 In the mathematical literature on open quantum systems, it is common to work from the beginning with baths in the 
thermodynamical limit; the bath state and dynamics are defined with the help of (a representation of) an abstract CCR algebra 
(see e.g. [20]). We shall not use such an approach here and treat together the infinite lattice limit for the particle and the 
thermodynamical limit for the bath. 
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An important property of our model is invariance under space translations. These translations are repre- 
sented on Tip by unitary operators Up(x) defined by Up(x)\y) ® |</>) = \x + y) ® \4>) for any x,y <E A and 
|0) € C N . We state some conditions on H^ op , which are in particular satisfied by the Hamiltonian (fj|. 

(Al) The hopping Hamiltonian H^ op has the form 

Ht P = Y, U £(- x ^°v U £& = E \y-^){z-x\{y\h hop \z) (5) 

where hhop is a A-independent self-adjoint operator on £ 2 (Z d ) £g) C N satisfying 

(y\hhop\z) = whenever \y\ > R or \z\ > R (6) 

for some R < oo. In particular, this implies that H^ op is translation-invariant, Up(—x)H^ op Up(x) — H^ op 
for any x £ A, and has a finite range independent of A. 

In the case of the discrete Laplacian ((H), we have hh op = (1/2) 2| z | = i(2|0)(0| — |0)( z l ~ l z )(0|) ® lc N - 

By using the Combes- Thomas estimate, which can be applied independently of S since the latter operator 
does not act on the translation degrees of freedom, one can show that the finite range condition (J6j) implies the 
propagation bound 

ll^le-^ly)!) < e ^ a W e -^-y\A ( 7 ) 

for some positive and A-independent constant k. Here and in what follows, ||A|| denotes the operator norm of 
the operator A (acting either on C^, Tip, or another space); in ([7]), the quantity inside the norm is a N x N 
matrix acting on the internal degrees of freedom of the particle. 



2.2 The bosonic baths 

The particle is coupled to one or several bosonic baths labelled by i € I (I is a finite set). Let T d be the 
d-dimensional torus, identified with the hypercube ] — tt, 7r] d . Let A* = (7r/L)Z d n T d \ {0} be the dual of 
the lattice A after having removed the origin^. The frequency Vi(q) of a boson with a (nonzero) quantized 
momentum q 6 A* is the value at q of a function i>i : q € T d H> Vi(q) G K+ (dispersion relation of the bath i). 
We assume that Vi is continuous on T d , C°° on T d \ {0}, and it satisfies 

Vi(q) >0 for q^O. (8) 

The Hilbert space of bath i is the symmetric bosonic Fock space built on f)i = £ 2 (A*), 

= F s (f) 4 ) = c e f)i e <8> s hi) e • •• (9) 

where <E) S stands for the symmetrized tensor product (see [28] for details and background). The full bath space 
is then given by Hp = ^> ieI V.p i . The boson Hamiltonian H A , acting on Hp, is 

ffB = EE Mq)< q a hq (10) 

where a* and ai iQ are the creation and annihilation operators for bosons with momentum q in the bath i. We 
recall that a* q and a^ q are unbounded operators on Hp, acting trivially on Hp j with j ^ i, which satisfy the 
canonical commutation relations [di, q ,a* ql ] — SijS qiq ' and [a,i,q,(lj,q'] = 0. 

2 In our finite- volume setup, bosons with zero quasi-momentum q = do not play any role and should be removed in order that 
all expressions in this section be well-defined. 
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2.3 Coupling between the particle and baths 

The particle and baths are coupled via a translation-invariant interaction Hamiltonian acting on *H A = Hp^Hp, 

= -|= E 9o,MWi ® e i? - x ® (a?, + a*,.,) (11) 

where X is the position operator (acting on £ 2 (A) as a multiplication by ir and acting trivially on C N ), Wi is 
an Hermitian N x iV-matrix that models the interaction with the internal degree of freedom, and go,i{o) are 
momentum-dependent coupling constants. Having in mind the thermodynamical limit which will be considered 
below, we assume that 30,1(9) are the values at the quantized momenta q G A* of some continuous functions 
go.i ■ T d H> C. These functions are called the form factors in the sequel. They must satisfy 50,1(9) = 9o,i(—q), 
q G T d , in order that i?A t be self-adjoint. Introducing the field operators 

<tf(vO = lAl" 1 / 2 J2 + PfaK«) (12) 



for (f G £ 2 (A*), one may rewrite if A as 



where we have set 

g x M = e [q ' x 9oAq) ■ (w) 

Up to the freedom in the form factors 30,1(9): the choice of the Hamiltonian (lll[) is dictated by the requirement 
that it must be invariant under space translations and linear in the creation and annihilation operators of the 
bosons. Space translations are represented on the bosonic Fock space of the bath i by unitary operators U A (x) 
satisfying U A (—x)ai^ q U A (x) = e lq x ai tq . One easily checks that for any x G Z d , 

U\-x)Ht,U\x) = Ht t (15) 

with U A (x) = Up(x)®(Q ieI U A (x). For instance, electrons in solids are coupled to low-energy acoustic phonons 
via an Hamiltonian of the form pT|) , see [5^1 [30] • 

The total Hamiltonian of the coupled system, acting on T~L A , is 

H t A ot = (X a Ht p + S)® 1 w a + l n ®H A + XH A t (16) 



where we have introduced the dimensionless coupling constant A in front of H A t . Using Vi(q) > and the 

J A 

z tot 



finitcness of A, one can apply the Kato-Rellich theorem to conclude that H A ot is self-adjoint on the domain of 



2.4 Initial state 

We assume that the particle and bosons are initially in a product state pp (£> pp, where pp is the initial density 
matrix of the particle and pp the initial density matrix of the bosons (i.e., pp and pp are positive operators on 
Hp and Tig with trace one). We now specify our assumptions on the boson initial density matrix pp. To this 
end, we consider the n-point correlation functions 

where af stands for a iyq or a* q . 

(Bl) The bath density matrix pp is translation-invariant and stationary with respect to the free dynamics 
generated by Hp : 

U A {x) P pU A {-x) = e-' ltH Bp A Q HH B = p A j Qr any x £ A, i £ I, and t G M. (18) 
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(B2) pg is quasi-free. That is, the correlation functions (|17p exist for any q±, ■ ■ ■ ,q n G A*, they vanish if the 
number of creators is distinct from the number of annihilators □ (in particular, if n is odd) and they satisfy 
the following Gaussian property (also known as Wick's identity): for n even, 

n/2 

HPB<] qi 4:, q ,---4:,j E n Hpb^At.J ■ ( 19 ) 

parings 7r of (!.,••• ,n) m—1 

The sum in U9\) runs over all pairing of (1, • • ■ ,n), i.e., all sets 7r = {(ti,t7i), • • • , (i„/2, c n / 2 )} of n/2 
pairs of distinct indices such that 1 = i\ < i% < ■ ■ ■ < L n / 2 , l"m < °~m f or an U m = lj ' ' ' i n /2, and 
{i>U"- , t n/2}U{cri,--- ,ct„/ 2 } ={1,--- ,n}. 

(B3) = ®i£ip% i is a product of quasi-free density matrices pg i on "H^ i . In particular, this implies that 
the two-point correlation function ti(p^ a f 9 1 1 a j'q 2 ) vanishes if i 7^ j. 

Note that, according to assumptions (Bl) and (B2), tr(pg a f q \ a fl 2 ) — also vanishes if qi 7^ q 2 . Assumption 
(B3) means that the baths are not correlated initially. Assumptions (B1-B3) imply that pg is completely 
determined by the set {Ci(l)',Q G A*} of occupation numbers (^(q) = tr(p^ a* „ai <q ) G R+ of bosons with 
momentum q in bath i. In particular, 

tr(^$fM$ffe)) = i^ (c A (g)^(g)^) + (i + C A (?))^)^ 2 (<7)) (20) 

for any cpi, ip 2 G £ 2 (A*). When taking the thermodynamic limit we will need the additional hypothesis: 

(B4) Ctio) are the values at the quantized momenta q G A* of some continuous function Q : T d \ {0} — > R+ 
such that |s ,i| 2 Ci € L 1 (T d ). 

The prime example of an initial state satisfying (B1-B3) is p^ = (£) ieI Pp. i where 

" S^K) (21) 

is a Gibbs state at inverse temperature /3j > 0. This situation corresponds to a particle coupled to thermal 
baths (which may have different temperatures). Then Ci(<?) = (e 13 ^ 1 ^ — is the Bose-Einstein distribution. 
If the form factors are such that Iffo.il 2 /^ € L 1 (T d ) then the last assumption (B4) holds true. 

2.5 Thermodynamic limit 

To observe irreversible phenomena, we have to consider the baths at the thermodynamic limit, that is, send 
A /■ 1 d keeping the boson densities fixed. For the particle, sending A /• 7L d corresponds to taking the infinite 
volume limit. 

Let i/hop be the bounded self-adjoint operator on "Hp = £ 2 (Z d ) (g) defined by (x\H^ op \y) — > (x\Hh op \y) as 
A /• Z d for any x, y G Z rf . This Hamiltonian is given formally by H\ wp — J2xei, d Up{—x)h\ lop Up{x). It describes 
the hopping of the particle between the lattice sites in the infinite volume limit. We identify all operators on 
£ 2 (A) (g) as finite-rank operators on £ 2 (Z d ) ® C N . By the finite range condition (Al), (x|i?hop|y) = for 
\x — y\ > 2R and H^ Qp —> Hh op strongly as A / Z d . Since H^ op (and thus Hh op ) are bounded, it then follows, 
e.g. by the Duhamel formula, that 

e -iw£ e -itH P strongly as A / Z d (22) 

with Hp = X a H hop + S. 

3 This property is sometimes called gauge invariance since it follows from the invariance of the state under the gauge transfor- 
mation a^q — > e di^q. 
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In what follows, we will denote by B(Wp) (respectively B%(Hp)) the Banach space of bounded (trace-class) 
operators on Hp (recall that the norm on Bi(Wp) is the trace norm ||A||i = tr(|A|)), and by Sp the convex 
cone of density matrices on B(Hp) (i.e., positive operators in B\(Hp) with trace one). We must clearly assume 
that the finite volume initial state of the particle converges as L — > oo in the trace-norm topology. 

(A2) p P -> PP in BiCHp). 

It is easy to show from the commutation relations of the af q 's that the field operator (|12p evolves under the 
Hamiltonian (flOl) as 

e ltH B^(ip)e' ltH B = <f>^(e u »*(p) . (23) 
The following space-and-time bath correlation functions 

/A(x, y; t, s) = tr(^$f (e^g x ^ (e isUi 9 v ,i)) (24) 

will play an important role in what follows. By translation invariance and stationarity of the bath initial state 
Pp (assumption (Bl)), f^(x,y;t,s) — f^(x — y,t—s) only depends on the position difference x — y and time 
difference t—s, where, according to (|20l) . 

//W) 4^ \9oM\ 2 Ut(Q> iq - x ^ tViiq) + (1 + CHu)>- iq - x e- u ^A . (25) 
By assumption (B4), f^(x,t) converges as A /* Z d to 

Mx, t)=J ^jL \go,M)\ 2 (CiWe^'e**™ + (1 + C^K^e"^) (26) 
uniformly in t (recall that v is bounded). 

2.6 Reduced density matrix of the particle 

The reduced density matrix pp(t) of the particle at time t > is the partial trace over Tig of the time-evolved 
density matrix of the total "particle + bosons" system, 

p P (t) =tr B (e- i *<tpA®p|e lt <t) . ( 27 ) 

The following proposition states that it is well defined in the thermodynamic limit under the assumptions 
described in the preceding subsections. 

Proposition 1 Assume that (A1-A2) and (B1-B4) are satisfied. Then for each t > and A > 0, the reduced 
density matrix (|27l) converges as L — > oo, 

p P {t) -» Z tA (p P ) (28) 

Ayz d 

in the trace-norm topology, where Z t ,\ is a completely positive trace- preserving map acting on Bi(Hp). 
This proposition will be proven in section 14.61 

3 Results and discussion 
3.1 The scaling limit 

To obtain rigorously a kinetic equation for ppit) = Z t \(pp), we perform a van Hove limit by setting t — A~ 2 r 
and letting A — ¥ while keeping the rescaled time r > fixed. In the interaction picture with respect to the 
internal Hamiltonian S, the reduced density matrix of the particle is in the scaling limit 

ftu» = lim e itS Z t ^(p P )e~ its (29) 

A-s.0,t=A- 2 T->co 
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(since we never consider objects on the bath space in this limit we write p s \(t) instead of pp, s i(t))- Note that 
the infinite volume limit (1251 has been taken first, before letting A — !• 0. Our main result is the existence and 
characterization of the limit (|29[) . Recall that A appears both in front of the interaction Hamiltonian H lnt and 
of the hopping term Hh op in the total Hamiltonian (|T6| . Hence hopping between the lattice sites goes to zero 
as A — > and the motion induced by the Hamiltonian -ffhop becomes effective only at large times t w A~ q t. 
It is then intuitively clear that for a > 2 the hopping will be absent in our scaling limit, whereas for a = 2 
both hopping and dissipative effects due to boson absorptions and emissions should be contained in the kinetic 
equation for p s i{r). 

Before stating the result, let us introduce some notation. In the following, {A, B} = AB + BA denotes 
the anticommutator of two operators A and B on Hp, {\s);s — l,--- , N} is the orthonormal basis of C N 
diagonalizing the internal Hamiltonian S, and E s £ cr(S) are the eigenvalues of S, that is, S\s) = E s \s). For 
any Bohr frequency u> £ cr([S, •]) = <r(S) — o-(S), we define the N x N matrix 

Wi,„= J2 (3°) 

s,s'=l,— ,N 

and the spectrally averaged hopping Hamiltonian 

< P = E S Es , Es ,( S \H hop \s%s)(s'\ (31) 

s,s'=l,— : N 

where S a ,b is the Kronecker delta symbol (equal to unity if a = b and zero otherwise). Finally, let us recall 
that a quantum dynamical semigroup (QDS) on Si('Hp) is a semigroup (7~ t )t>o of completely positive maps 
T T : Bi(Hp) — > Bi(Hp) preserving the trace and such that t 6 M + h-> T T is *-weakly continuous. Lindblad [T] 
has derived the general form of the generators of norm-continuous QDS (see also [2] and an extension to 
unbounded generators in [3"T]). 

Theorem 1 Let assumptions (A1-A2) and (B1-B4-) be satisfied and let a > 2. Assume moreover that the 
infinite-volume correlation functions fi(x,t) satisfy 

1 f°° i*i 
lim - / dt sup \fi(x,t)\e- » = . (32) 

Then for any r > and any pp € Sp, the limit (|29j) exists in the trace-norm topology and is equal to p s i(r) = 
e rC pp, where (e rC ) r >o is a norm- continuous quantum dynamical semigroup on B\{J-Lp). If a — 2 the generator 
& of this semigroup is given by 

C\p) = -i[< p + T, P ]+ E( E City-x,")Wi, u ®\x){x\pW* u ®\y){y\ 

-^j^{w* u Wi,u l t * m , p) J (33) 



where 



T= E 3 {/ dt /;((), t)e- iiw j WlJVi,*®!^ (34) 



uea{[-,S]) iei 

and Ci(x,ui) is the time Fourier transform of fi(x,t) 



a{x,w) = / dtf^x,^-^ . (35) 
Ifa>2, is given by the same expression as in Ii33\) but without the term — i^hopj/ 3 ] - 
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To see that £)■ in ([33| has the Lindblad form (and therefore that (e rC ) T >o is a QDS), we first rewrite 

^ (P) = -i[< P + T, p] + St(p) - \{9l*{l),p} (36) 

where the map 21 abbreviates the term involving a sum over x and y in the right-hand side of (1331) and 21* is 
its adjoint with respect to the trace, i.e., tr(2l(p)A) = ti (p9l* (A)). We note that Ci(x,uj) is of positive type in 
the x- variable (this follows from the fact that fi(x,t) is a correlation function, see (|24p) and therefore it is the 
Fourier transform of a positive measure Ci(dq, uS) on T d : 

c i (x,cj) = -^ [ ci(dq,u)e-^ x . (37) 

This shows that 21 has the Kraus form 




with Vi^ = (2n) d / 2 ®e lq ' x . Thus 21 is a completely positive map [32 . Moreover, 21 is bounded on B\(T-Lp) 
because for any p £ B{Hp), p > 0, 

||a(p)||i = tr(a( P ))< E^(°' w )II^H 2 IIHIi (39) 

w£<r([-,S]) ie/ 

and Ci(0,w) is finite by the integrability of the correlation function fi(0,t) (assumption (1521 ). Since also 21* (1) 
and i?ho P + T are bounded operators on Hp, it follows that £ 1 ' is bounded on Bi(Hp). This boundedness and 
the completely positivity of 21 imply that the operator in (|36|) generates a norm-continuous QDS [1]. 

Another way of phrasing Theorem [T] is to say that the rescaled density matrix p s i (r) satisfies the Bloch- 
Boltzmann master equation 

±p sl ( T ) = £\p sl (r)) . 

Note that commutes with i[S, •], a generic fact for generators obtained via a weak coupling limit [3 [3]. The 
self-adjoint operator T in acts trivially on t 2 (X d ) and commutes with S; it represents the energy shifts of 
the particle due to its coupling with the bosons (Lamb shifts). In the following sections 13.21 and HOI we unwrap 
the form of the generator £} in two different situations and discuss the physical phenomena described by the 
corresponding master equation. 

The major technical assumption of Theorem [1] is the integrability condition f|32[) on the boson correlation 
functions. This assumption should be compared to the analogous condition for confined systems [31[7J[H], i.e., 
the integrability of the correlation function (|26p for x — 0. An explicit computation and the use of a stationary 
phase argument performed in Appendix IA1 yields: 

Proposition 2 Let us assume that the form factor go^ has a support contained in the open ball {q £ T d ; \q\ < 
7r}, that |go,i|(<z) an ^ Ci(o) depend only on the modulus \q\ of q, and that the bosons of bath i have a linear 
dispersion relation Vi{q) — \q\ on the support o/<7o,i- Furthermore, let the functions if)i^(\q\) = \go,i(q)\ 2 Ci(l) 
and ipi-(\q\) = \g ,i(q)\ 2 (l + d(q)) belong to C 2 (]0, tt]) and 

M min{d - 3 '^V±(M) , M'-VUG?!) , M d "X±(M) ( 40 ) 

be integrable on [0,7r]. Then assumption (3ty) on the correlation function fi{x,t) is satisfied in dimension d>2. 

If the bosons are initially at thermal equilibrium, so that Q(q) = (e' 3 ' 9 ' — the assumptions on ipi,± m 

Proposition [5] are satisfied if go,i{\l\) e C 2 (]0: 7r ]) an d 

|?|min{d -4,^ }| |2 | |*-8 d a i&nd rf _ 2 2 

d\q\ d\q\ z 

are integrable on [0,tt]. 
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Remark 1 For most natural models, assumption ¥3ty fails in dimension d = 1; see AvvendixXM for an explicit 
computation. 

Remark 2 For physical applications it would be of interest to estimate the error terms in the convergence to 
the scaling limit, that is, to have an explicit bound on \\pp(\~ 2 r) — e~ lA t ^ S ''^p s i(t)\\i. Due to the repeated use 
of the dominated convergence theorem, our proof of Theorem [7] does not exhibit such a bound. One could in 
principle obtain the error terms by assuming some explicit decay of the correlation functions fi(x,t) as in fffffl 
(see Appendix B in this reference). 

Remark 3 The choice of the scaling exponent a in the factor \ a in front of the hopping Hamiltonian i/hop 
in Mb)) is more dictated by mathematical than by physical motivations. For a > 2 the hopping Hamiltonian 
-ffhop is absent in the dynamics in the scaling limit and one has a trivial coupling between the hopping and the 
internal degrees of freedom in the Lindbladian For a < 2, as explained in the introduction, the convergence 
in the weak coupling limit A — > is mathematically more involved. It is not clear whether in this case one can 
still distill a limiting Lindblad operator. 

As already indicated in the introduction, in the case a — 2 our scaling limit incorporates both features 
of the singular and weak coupling limits: the translational degrees of freedom are treated within the singular 
coupling limit and the internal degree of freedom is treated within the weak coupling limit. The fact that the 
Hamiltonian 

H x = X 2 H hop + H B + XH int (42) 

corresponds in the limit A — > to the singular coupling limit was pointed out by Palmer the dynamics 
generated by (|42l) at time t can be mapped into the dynamics generated by the Hamiltonian 

H' x = H hop + H' B + Hfa t (43) 

at the (unrescaled) time t and with an (unrescaled) hopping term Hi lopi where H' B is a free boson Hamiltonian 
as in (jlOp . The new interaction Hamiltonian H^J is not multiplied by A and is given by the original interaction 
Hamiltonian -ff int , as given by (|13[) . but with a rescaled form factor g\2 x (X 2 q) instead of g x {q). The two dynamics 
generated by (|4"2"|) and (|4"3")) are exactly the same, in the sense that e ltHx and e 1THx coincides up to a conjugation 
by the unitary operator transforming the annihilation operator a(<p) of a boson with wavefunction ip e L 2 (T d ) in 
the initial problem onto the annihilation operator a'(cp\) of a boson in the new problem, with ip\(q) = Xip{X 2 q). 
In the limit A — > 0, the rescaled form factor becomes singular. 

Let us mention that a model of a quantum system coupled to a free fermion bath has been studied in |23j 
in the singular coupling limit, using the approach of Refs. [7J [5]. 



3.2 A jump process in momentum space 

Let us choose -f/hop = — A according to ((4]), a — 2, and assume that all baths are initially in thermal equilibrium 
with the same temperature = fi^ 1 ■ We may then just as well consider only one bath initially in the Gibbs 
state p B = p^. 

This setup is interesting if one sets out to study diffusion and decoherence of the Brownian particle, see 21j. 
The resulting master equation can be written as 



= -i[A + T,p sl (r)]+ £ J A c(q,u)(w u ®T qPsl (T)W:®T; 

--{w^Sl^.puCr)}) (44) 

where T q = e lq X is the unitary momentum translation operator on l 2 (jL d ) and we have used the notation 
Ci(dq, lu) = c(q, uj)dq with c(q, lu) the positive distribution which for lu ^ is given by (see (|2T))) , (l3"5j) , and (|37l0 

^q,u) = 2n\go(q)\ 2 (((-q)5(v(-q)-u J ) + (l + aq))5(v(q)+Lo)) . (45) 
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Here £(g) = (e^"«' — is the Bose-Einstein distribution introduced in (f2T)|). and we have assumed that 
Ci(dq,uj) defines a bona fide measure, which requires some mild additional conditions on v. The first term in 
the rate (|45[) corresponds to absorption processes of a boson with momentum — q and frequency v(—q) = ui] 
it is proportional to the mean number C(~<z) °f bosons with momentum — q. The second term corresponds to 
spontaneous and stimulated emission of a boson with momentum q and frequency v{q) — — w and is proportional 
to 1 + (the term 1 comes from spontaneous emission). The delta distributions in (1451) accounts for energy 
conservation (see Figure [T]). Note that we have the detailed balance condition c(q,ui) = c _ ' 9lj c(— q, — oj). 

To appreciate the equation (|4"4")l , it is worthwhile to see what it implies for the evolution of diagonal elements 
of the density matrix in the eigenbasis of S and the momentum basis for the translational degrees of freedom. 
Let us define (assuming that the sum on the right-hand side is absolutely convergent) 

p sl (r;k,s)= e-^- x '> k {x, S \p sl (T)\x',s) (46) 

with \x,s) = \x) ® \s). The momentum density (j46|) satisfies J2^=i It* dfc Ps\(t; k, s)/(2Tr) d = 1 for any r > 0. 
For simplicity, let us assume that the spectrum of S is non-degenerate, i.e., all eigenvalues of S are simple. Then 
(gH) gives 

f) N r 

— Psl (r; k', s') = J Td dk(-f(k\ s'\k, s) p sl (r; k, s) - 7 (fc, s\k', s') p sl (r; k', s')) (47) 

where 

-f(k',s'\k,s) = {2n)- d c{k' -k,E s , - E s )\(s'\W\s)\ 2 . (48) 

In formula (|47[) one recognizes the structure of a forward Markov generator with (singular) transition rates 
7(fc', s'\k, s), acting on densities of absolutely continuous probability measures (hence on L 1 -functions) on T d x 
{1, • • • , N}. Therefore, the master equation (l4"4l describes the stochastic evolution of a particle with momentum 
k and internal state s, which may jump from the state (fc, s) to (k',s') by emitting or absorbing a boson of 
momentum q and energy v(q), as represented in Figure [T] According to (|45[) and (|48l) . only jumps satisfying 
energy and momentum conservation v{q) = \E s i — E s \ and q = sign(£v — E s )(k' — k) are allowed (here sign is 
the sign function). Note that, in the limit A — > 0, the energy of the particle coincides with its internal energy 
E s since the hopping energy was assumed to be of the order of A a with a > 2. This explains why the detailed 
balance condition 

7 (fc', s'\k, s) = e^ E '- E ^j(k, s\k', s') (49) 

does not involve the kinetic energy but only the internal energy levels E s and E s i . 

In the absence of internal Hamiltonian (i.e., for 5 = 0), energy and momentum conservation and the fact 
that v(q) ^ for g / imply that, up to second order in A, the particle can only emit or absorb zero- 
momentum bosons and thus cannot change its momentum. Hence the coupling with the bath has no effect on 
the translational degrees of freedom in the scaling limit ([29)) . This features would not change if we consider 
baths at different temperatures, or a more complicated hopping Hamiltonians. 

Given some mild technical conditions, one can show that a particle described by the Lindblad equation (|44"]) 
diffuses in space, i.e., J2x£Z d \ x \ 2 ( x \Ps\( T )\ x ) oc t as r — ^ oo, see [2lT 122] . 

3.3 Ratchet 

For this application, we choose 

ffhop = -^2 (\ x )( x + e i\ ® I^X-M + |z + ei)(x| <E> |»«-><fi_>|) • (50) 

Here, \s^) and |s<_) are two distinguished eigenstates of S and we singled out a spatial direction by coupling 
these states to the motion in the direction of the unit vector e\ £ J/ 1 . We assume that S has four distinct 
eigenstates labelled by s-j-, sj_, s_,., and satisfying 

S|« t ) = |s t ), S|s;)=-| S; ) , and S|s_>) = S|*<_) = . (51) 
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(a) The particle makes a transition (k, s) — > (k' , s') and (b) The particle makes a transition (k, s) —> (k 1 , s') and 
emits a boson of momentum q and energy v(q) with absorbs a boson of momentum q with k + q = k' and 
k = k' + q and E s = v(q) + E s , . E s + u(q) = E 3 , . 

Figure 1: The processes contributing to the gain term (first term on the right-hand side of H47I) - ), Emission corresponds to 
E s > E a i and absorption to E s < E a i . 

For simplicity, we choose equal eigenvalues of S corresponding to the states |s->) and |s<-), so that H^ op and 
S commute. We will exploit the fact that we have two reservoirs. Bosons from the first reservoir couple to the 
s-variable by 

W i=1 = | S t)^| + |s 4 )(^| + |0(stl + M^| (52) 
while bosons of the second reservoir couple as 

W %=2 = | St )^| + + |*<_)fol + l s ->WI • ( 53 ) 

We choose the initial state of the reservoirs to be pb — p^ <8> , where ui^ is a Gibbs (thermal) state at 
temperature T; = f3~ l ■ If the two reservoirs have the same temperature, then the model does not display any 
current. However, by preparing the reservoirs at different temperatures T\ ^ T 2 one breaks the time-reversal 
symmetry and a current will in general emerge. To simplify the forthcoming discussion, we choose go.i = ffo,2 
(the form factors are equal) and 

Ti = 0, T 2 > . (54) 

The boson field induces jumps between the internal states \s) as represented in Figure [2] By using energy 
conservation (see (|45jl ). one easily convinces oneself that, with the temperatures chosen as above, the particle 
can make a transition from |s<_) to |s->) by emitting a boson of the first reservoir and absorbing one of the 
second reservoir, whereas it can not make the reverse transition from |s->) to |s«_) (there are no boson with 
frequency to = 1 in the first reservoir at T\ = 0). Since all the jumps between eigenstates of S happen at 
fixed position, these transitions do not in themselves induce a current. However, the hopping Hamiltonian 
_ffhop allows for transitions between the states |s_>) and |s<_). Hence, a current flows in the ei-direction. The 
possibility of extracting work from the system is already visible in Figure [21 where one sees that the particle can 
go from |s<_) to |s_>) via the upper (respectively lower) level only clockwise (anticlockwise). Once one has this 
property (which is excluded in equilibrium), it is clear that one can devise a scheme to convert this "internal 
current" into a spatial current. 

4 Proofs 

4.1 Preliminaries 

Our proof of Theorem Q] follows a similar approach as in some previous works of one of the authors, in particu- 
lar [2T1 [34] . The starting point is a Dyson expansion of the propagator Z ti \. We consider the situation in which 
the particle is coupled to a single bath. Therefore, we do not write the lower index i and the corresponding 
sums, which do not play any role and obscure the notation. The proof for several baths, i.e., |/| > 1, is exactly 
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k> 

Figure 2: Possible jumps between internal states of the particle. Note that the jumps due to the coupling with the reservoir at 
temperature T\ = can take place only from higher to lower energies, whereas the jumps due to coupling with the other reservoir 
at positive temperature can go either way. 

the same. We also restrict ourselves to the case of a scaling exponent a = 2. The case a > 2 is simpler and can 
be treated similarly. 

Let us first fix some notation. For products of operators, we use the conventions 

<- -> 

Y[ Aj = A n ... A 2 Ai, Yl A J = A i^2 ...An. (55) 

j — l....,n j— l,...,n 

The operators acting on Bi(Hp) are denoted by calligraphic fonts, e.g. T or £, and we use their (operator) 
norms 

HT||= sup Mi. (56) 

AeBx{Up) \\ A \\l 

We denote by C constants that can depend only on the space dimension d, the internal space dimension N, 
the parameter k in the propagation bound (J7J, and the correlation function f(x,t). 



4.2 The Dyson series 

7A 

7 A 



Let be the propagator for the reduced dynamics of the particle at finite volume, defined by 



04) = tv B (e-^- (p P ® p%) e*<« ) , p P eS(H A P ). (57) 



We will show below that the Dyson expansion with respect to the interaction Hamiltonian ((131) °f this propagator 
converges in norm; this expansion reads 



n>l J0<ti<~<tn<t (x u - ,x n )EA» 

ti B ( [v A (t n ) ® (t„) , • • • [v A (t 2 ) ® < (i 2 ) , (*0 ® <(ti) , ^ ® ^]] • • •] ) (58) 
where <I> A (i) = $ A (e itiy <7 X ) is the freely-evolved field operator, see (|23[) . and 

\/ A (t) = e wi ^W® |x)(a ; |e- itff i (59) 
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For any xeZ d ,t>0, and T G B(%j,), let us set 

. V x A (t)T \£l = L (f A (x,t) i(l = L 

2 A (x,t,l)(T) = { , f A (x,t,l) = { (60) 

-TV x A (t) if I = it I / A (^,t) if/ = i? 



where f (x,t) is the correlation function (|25|) . Note that f A (x,t) = f A (—x,—t). By using the quasi-freeness 
assumption fS^, one gets 

^ = 1 +E E E / d*^„(7r,i,x,i) (6i) 

n>l pairings 7r (x,i)£ A 2 " x{L,i?,} 2 " 0<ti<--'<t2 n <t 

where di stands for dti ■ ■ ■ dt 2 m the sum runs over all pairings n — {(ti, cri), ■ ■ ■ , (i n , c„)} of (1, • • ■ , 2n), and 

<— n 

V^(7L,|,^I) - (-A 2 )" J] I A {xj,tj,lj) J] / A (x CTm - z tm ,i CTm - t Lm ,h m ) (62) 

j — 1,--- ,2n m— 1 

if xi, ■ ■ ■ ,X2n € A, and V A n (n, t, x, V) = otherwise. We do not write explicitly the dependence of Hp, V A , 
and X A on the coupling constant A to simplify notation, but we keep it in V A n and D A \ because we will later 
consider the limit A — ¥ of these quantities. 

Already at this point, we can establish the norm-convergence of the series ([BT]). Indeed, \\l A {x,t,l)\\ < \\W\\ 
and \ f A (x, t,l)\ < / A (0, 0), thus the nth term in the series (|6"Tj) has a norm bounded by 



t 



2)i 



(2|A|A||^||) 2 "^)y(/ A (0,0)r 



((|A|A||^p) 2 2/ A (0,0))" 



(63) 



where the term between the square brackets [•] is the number of pairings 7r of (1, • • • , 2n). Hence the Dyson 
series (|6ip at finite volume converges in norm. One can prove that its sum is equal to the propagator in the 
interaction picture, 

V A X = jWUz A x . (64) 

We will consider the Dyson series at infinite volume and we simply drop the superscript A on Hp, I, V, and 
/ to denote the corresponding objects for A = Z d . We argue that in the infinite volume limit A f 

e -tt[»£.-] _> e -it[Hp,-] } l A (x,t,l)^l(x,t,l), V^ n {TL,t,x,r)^V x ,n{K,t,x,l) (65) 

strongly on B\(Hp). Indeed, recall that given some bounded operators A A on a Hilbert space such that A A — > 
and (,4 A )* -> strongly, then P A T||i and ||TA A ||i for any T G B\(Hp). The first limit in ([65) 
follows from this property, the inequality ||j4 a TC/||i < ||v4 A T||i for U unitary, and the strong convergence of 
e - itH p on Hp, see 1(22]). As V A (t) V x (t) strongly on H P (again by (|22|))). the second limit follows from the 
same property. Since V A n (n, t, x, I) is a finite product of the X-operators and correlation functions / A , the third 
limit in (|6"5"|) then follows from the second one and from the pointwise convergence f A / (see Section dS]). 
This implies that term-by-term, the series "D A X converges strongly to the infinite-volume Dyson expansion 

^^EEE / dt - V *AlL,l2L,L) ■ (66) 

™ 2L x,l Z2n(t) 

The convergence of this series has not been addressed yet, so for the moment we consider it as a formal series. To 
prove Proposition [TJ we will show that this convergence is in some way uniform in A and apply the dominated 
convergence theorem. We have made in (|66[) the following abbreviations, which will also be in place in the 
remaining of the paper: 

• we sum over n = 0, 1, . . ., where it is understood that the term corresponding to n = is equal to 1; 

• the sum over 7r ranges over all pairings of (1, 2, • • • , 2n); 

• Z n (t) denotes the simplex {t = (ti, ■ ■ ■ , t n ) £ [0, t] n ; t\ < t2 < ■ • • < t n }; 

• the sum over x and I range over Z 2nd and {L, R} 2n , respectively. 
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(a) 



(b) 
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Figure 3: (a) crossed diagram (b) ladder diagram. 
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4.2.1 Graphical representation 

It is convenient to represent Va,ti or V^ n by a diagram in which all times t\ < ti < • • • < ^2n are ordered on 
the real line and pairings are represented by bridges linking two distinct times. Two examples of diagrams are 
represented in Figure [3] Replacing the times by their indices 1, 2, • • • , 2n, the diagrams with 2n points are in 
one-to-one correspondence with the pairings of (1, 2, • • • , 2n). A diagram containing two pairs (i, a) g 7r and 
(t',cr') g 7r such that i < i! and l' < a is called a crossing diagrarr^. A non-crossing diagram will be called a 
ladder diagram; it corresponds to the pairing 7r laddor = {(1, 2), (3, 4), • • • , {In — 1, 2n)} respecting the order, see 
Figure El 

The strategy of the proof of Theorem [JJ consists in showing that in the scaling limit A — > 0, t = A~ 2 t — > oo, 

• all crossing diagrams in the infinite-volume Dyson series (|66j) converge to zero; 

• the ladder diagram in (|66j) of order 2n converges to the corresponding diagram of the Dyson expansion of 



'e rC , where & is the Lindblad generator given in ([33 



4.3 Plan of the proof 

Below, we give the main steps of the proof of Theorem [TJ 
4.3.1 Topology 

We first introduce a notion of convergence that is particularly useful for the problem. Let (7a) a be a family of 
operators on Si (Hp). We associate to 7a a kernel with values in operators on the vector space B(C N ) of N x N 
matrices, defined as follows: 



(Tx) Xo , yo ; X , y (M) = {x\Tx(M®\x )(y \)\y) , M g B(C N ) . 
We write pt-lim A ^ 7a = whenever 

lim ^2 \\(7x)x ,yo;x,v\\ = for any Xo,yo g Z d , 



(67) 



where the norm inside the sum is the matrix norm. 



Lemma 1 Let 7a be uniformly bounded operators on Bi(Hp). If pt-lim^g 7a = then 7a — > strongly, that 
is, limA^o ||7a(T)||i -> for any T g Bx{% P ). 



Proof. Note first that for any T g Bi(%p) 

MTA)\ 



\ T \\i = sup 

AeB{H P ) \\ A \\ 



< sup £ tr CN (\( X \T\y)\)M^R< N £ Ux \ T \ 
AeB(H P ) x ~^ 7J i \\ A W ^ 



(68) 



4 Traditionally, some of these diagrams are called "nested" , but we call all of them "crossing diagrams" as we represent them on 
a single time axis, see Figure \3\ 
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where the last inequality follows because tr(|M|) < A||M| for any hnite matrix M e B(C N ). Let T have 
finite support in the sense that (xq\T\uq) is nonzero only for a finite number of Xo,yo G 7L d . Then, by (|68l) . 
pt-lim A ^ 7a = implies || (^) 1 1 1 —> 0. Next, one checks that operators T with finite support are dense in 
Bi(Hp). Actually, let Pa = YlxeA \ x )( x \ ® lc s be the finite-rank projector on span{|a;);a; € A} ® C^, with 
A = 7L d /(2LZ) d as before. Without loss of generality, we may assume that T > 0. If is an orthonormal 

basis of "Hp diagonalizing T and Pj > are the eigenvalues of T, then 

||PaTFa - Til, < ^^(IIPaI^X^KPa - l)^ + || (Pa - Wi){1>j\\\i) 

3 

< 2£>||(P A - 1)|^}|| 

3 

as L — > oo by dominated convergence. Hence, for any T £ Bi(Hp) and e > 0, one can choose a decomposition 
T = T +T 1 such that T has finite support and ||Ti||i < s. Then ||7x(T)||i < ||7a(T )||i + Ce by the uniform 
boundedness of (7a) a • The claim follows. □ 

A consequence of lemma [1] is that in order to prove the strong convergence of T>\-i T \ to T> T in the scaling 
limit, it is enough to show that 

pt-lim {V X -2 T a - V T \ = . (69) 

A-s-0 

In fact, one has ||2?*,a|| = 1 for any t by the following standard argument. One first notes that the finite 
volume propagator Z^ x defined in (|5T[) preserves positivity and trace and hence so does the sum 

v t.\ of its 

Dyson expansion (related to Z t A A by (1M10 as well as its strong limit T>t.\ as A | Z d (we assume here that 
Proposition [T] has been already established). Then the dual of T> t ,\ under the trace, T>* t A , which acts on B(Hp), 
preserves positivity and satisfies X>* A (1) = 1. It follows that \\T>^ A || = 1 (see e.g. PS, Corollary 3.2.6) and thus 

PM = i. 

4.3.2 Assumptions 

In the remainder of this paper, we always assume (A1-A2) and (B1-B4) to be valid without further mentioning 
it (note that once the Dyson series for Z t ,\ is accepted as the basic object of study, one does no need those 
assumptions anymore). To prove Theorem [TJ we rely on: 

a) The propagation bound ([7]), but for A = Z d : by (|2"2"|) and as the right-hand side of is independent of 
A, this bound remains valid for A = Z d (alternatively, one can check this directly by the Combes-Thomas 
estimate, using the fact that -ffhop has a finite range). 

b) Assumption (I32|) on the infinite volume correlation function f{x,t) or, in most intermediate steps, the 
weaker requirement that f(x, ■) is integrable for any x E Z d . 

4.3.3 Step I 

We first prove that the Dyson series 2? A -2 T A , considered as a series in n, x, and I, converges absolutely and 
uniformly in A, in the sense that 

EEEE su p / dt\\(v x , n ( K ^,xd)) XO!yo .j\<™ (70) 

n x,I tt x,y X>0 Z J (X -, T) 

for any Xo,yo G Z d (see Proposition |3] below). This bound is the crucial point in the proof, since it allows us to 
estimate the perturbation series term by term. It relies heavily on the assumption (|32p . 

By a similar bound on the finite- volume Dyson series 2^\> the ^ ac ^ that this series converges term by term 
as A j" Z d , Lemma [T] and (|65|) . we obtain Proposition [TJ with 

Z t ,x = e- H ^V t , x . (71) 
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Here the integrability of / A (x, •) is not needed, the only requirement is its pointwise convergence as A | Z d 
(which follows from (B4))- 

These results are accomplished in subsections 14.41 and 14.51 

4.3.4 Step II 

We show that every single crossing diagram vanishes in the scaling limit, in the sense that, for any n > 1 and 

xo,yo,x, y e Z d , 

lim / dt\\Vx, n {TL,i,x,V)x o , Vo ;x,v\\=0 whenever n ^ 7r laddcr . (72) 

A— >0 J 

This step is essentially taken over from the original work [7] , but we review it in subsection 14.71 To feel why 
this holds true, note that when t = \~ 2 t 3> 1, the time integration domain of a crossing diagram of order In 
is much smaller than that of the ladder diagram of the same order. This is due to the restriction t b < t L > < to- 
associated to the nested pairs (i, a) and (i',a'), see Figure El 

By dominated convergence and (|70|) . this implies that the contribution of crossing diagrams in the Dyson 
series 2?\- 2 t a vanishes in the limit A — > in the topology introduced in section f4. 3. 11 i.e., 

EEE E / m^2L,t,x,l)) xoMx J ->0 (73) 



n v,l x,y 7L#7r ladderZ2n(A _ 2T) 



for any Xo,yo £ Z d . 



4.3.5 Step III 



It remains to evaluate the contribution of ladder diagrams 7r ladder . Let IC(t,x,V) be the bounded operators, 
defined in Proposition [Bj below, which yield the limiting QDS (e rC ) r >o upon summing over n, x, I, and r: 

E / dr^/C n (r,^) = e^ (74) 

" Z n (t) -'- 

where the sums and integrals are absolutely convergent in norm. We will show in sections 14.81 and 14.91 that 

lim / dt VA,n(7r ladder ,i,x,{) = / drlC n (T,x,l) (75) 

A->0 J J 

in norm on B(Bi(Hp)), which of course implies 

J dt^n^feaO^^^ / d L (JC n (T,x,l)) x(hya . xy . (76) 

Z 2n {\- 2 T) Z 2n (r) 

Dominated convergence allows us to conclude from ([701) and ([75]) that holds true. 

4.4 Estimating each term of the Dyson series 
Lemma 2 Fix n,n,xo,yo, and t€ Z2n(t). Then 

n 

\\(VxUK,t,x,L)) XOtyo . !X J < (CA) 2 "e 4 ^ J] h n (t am ~t L jR x ^ x Jx,l) (77) 

m— 1 

where 

h n = sup \f(x,-)\e-^ (78) 
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Figure 4: Equivalent representations ofan = 8-points diagram: on the right diagram, the times tj with lj = L 
are put on the upper axis and the times t k with l\. = R on the lower axis. Here |Z| = 5, (ji,j2i33i3A, 3s) = 
(l,3,4,6 ) 8),and(fci,fe 2 ,* 3 ) = (2,5 s 7). 



and Rx^]y ;x,y{x,V) is independent of X and A and such that 

n -° x,L x,y 

The bound J 77] ) is also true if one replaces V n by and h n by max xe A |/ A (x, •)!■ 



(79) 



Proof. We hrst give an explicit formula for Vx >n (w,t, x, f)(M (8 \xo)(yo\) in terms of the operators V x (t), with 
M £ B(C N ) and xo, yo G 7L d ■ For a fixed I £ {L, R} 2n , let |Z| denotes the number of indices j such that lj = L, 
j = 1, • • ■ , In. We change the labelling of the indices and coordinates by defining (see Figure @| 

{ji<---<J|i|} = {je{l,-.-,2n}; lj=L) 

{h < ••■ < fc 2 «-|i|} = {k 6 {!,••• ,2n} ; l k = R) 



and 



xo 



Xo 



Xl 



Xk± j 



' ^1 
i 2/2n-|J| 



•'-'A- 



an-m 



It follows from ([60| and J62j that for any matrix M e ^(C^), 

VA, n (7r,t I x,I)(M®|xo)<2/o|) = A aw (-l) w+ l | 'V^ |l| (*,-,„ ) • • • % fajM <8> |x >(y | 

n 

^(^•••V^-l^fcan-l.i) II ~*i»>0- (80) 

m— 1 

Note that a similar formula holds at finite volume for V^ n (n,t,x,l). 

We must bound the norm of the right-hand side of (|50)l . Let us denote by G(t;x,y) = (x\e~ ltHp \y) the 
time-dependent Green function associated to the free motion of the particle at infinite volume. Using (|59[) and 
setting tj a = tk =0, one gets 



'l^ii|(*iii|)---^ I (*j 1 )l a; o) = (?(-*,•„! ;x,X|i|) Y[ WG ( t j P -tj P -x\Xp,Xp-i) 

P =i,...,\i\ 



(81) 



and 



>y2n-\i\>y) ■ (82) 

g=l,...,2n-|/| 

Thanks to the propagation bound ([7]) , we have 

||(V A ,„(7r,t,x,/)) II (83) 

\l\ 2n-UI n 



< (||VF||A)^ n e' 



2n„4A^t, 



\x-Z\i\\-\y-ii2 n -\i\\ TJ c - 
p=i 



1 n e- 

9=1 



|y<j-yg 



- 1 ' II MW-Oe" 
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Next, observe that 



|2| 2n-|i| 



E I £ p ~ + E I&2 ~ y-J-il > - E \ x<r ^ ~ x <-m\- l^o - yol • (84) 



n 

p— 1 g— 1 m— 1 

Actually, (|84[) is a consequence of the inequality 

2n+l 

^ \z m +i - 2m | > max{|z CT - z t | ; t, a = 0, ■ ■ • , 2n + l} 

m=0 

applied to (z , • • • , Z2n+i) = (^|i|,--- ,xo,yo,--- ,V2n-\i\) € Z (2ll+2)d . Replacing (JM]) into (JHHJ) , one gets the 
result with 

-2n 2n-\i\ 
Rx„, yo ;x,y{x,l) = (2 ^ V |:E " £ ' i ' 1 ~ l ^ 2 "- 1 1 ' 1 e^ l*o-ifti I e - ' 1 1 JJ e"' 15 ' - "'- 11 . (85) 

zez d p=i 9=1 

The proof for the finite lattice is the same, since we only used the propagation estimate. □ 



4.5 Uniform convergence of the Dyson series 
Proposition 3 Let Xq,uq £ Z d . For fixed A and t, one has 



EEEE / d * 8up ||(V^ B (7r J t > x > i)) a!o ||<oo (86) 



n x,l K *,V Z2n{t) 



where the supremum is taken over A =] — L, L] d D Z d for all finite L > 0. Similarly, assume that the correlation 
function f(x, t) satisfies (|32|) and fix t > 0, then, 

n £l s,y > tt Z2n(A -2 T) 

Proof. We first show the second claim. The proof of the first one follows similar lines. We bound ([57)1 with 
the help of Lemma [2] by 

~ n 

Co, Wo e 4r E su p{( CA ) 2 " E / d ^ II M*»»-0} (88) 

" A> ° = Z a »(A-»r) m=1 

where Ck VQ < oo denotes the supremum in (j79[) . The sum over all pairings n and the time integrals are 
conveniently rewritten with the help of 

Lemma 3 For any locally integrable function g : R 2 ™ — > R and < to < t, 
E J d * 5 (*n ; ; • " - ;*l„)*<7 n ) 

— t <tl<-"<t2n<t 



y dni---dit„ y dui • • • du^sf(u l! u' 1 ; • • • ;u n ,u^) 
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We leave to the reader the proof of this lemma, which is based on change of variables. 

Applying Lemma[3]with g(u±, u[; ■ ■ ■ ; u n , u' n ) = Y[ m hn(u' m — u m ) as in (|55)) . bounding the integrals over the 
u m and u' m by ||/in||i — (fn° dth n (t)) n times the volume of the n-dimensional simplex Z„(A _2 r), we conclude 



that (|88p is smaller than 

a,o, w e 4r £^(<? 2 ||MiT) n • (89) 

n>0 

By Stirling formula n\ ~ (27rn) 1 ' 2 (n/e) ri as n — > oo, one finds that the convergence of the series in (|89|) is 
ensured by assumption (|32l) . that is, by the condition ||/i n ||i/n — > as n — > oo. Thus the second claim of 
Proposition [3] is proven. 

To show the first claim, we replace h n by max Ig A |/ A (a;, -)| (see Lemma[2]) in (|88|) . Here we do not need to 
assume that this function has a finite I^-norm, we bound it by / A (0,0), see (|25|) . The quantity in ([86)) is thus 
smaller than 

1 /C 2 A 2 i 2 sup A / A (0,0) N " 



^o,y e 2^ 3 o < 00 ' ( 90 ) 



n! 

Note that sup A / A (0, 0) is finite since / A (0, 0) converges as A f Z d . This concludes the proof of Proposition [3] 
□ 

4.6 Proof of Proposition [I] 

One has 

||((V A , n -VA,„)(2r,t,x,/)) :co!!/o . a;) J| — > as A / 7L d . (91) 

Actually, let us set A A = (V A n — Vx tn )(TT,t,x,l). Since {A A ) Xo>yo - Xiy is a finite matrix, it is sufficient to prove 
that \\{A A ) xo<yo - x , y {M)\\ < II (A A )x ,y ;x,y(M)\\ i for any M € ^(C^), so that the convergence dHJ) follows 
directly from the last claim in (165[) . By Proposition [3l (|9ip . and dominated convergence, for any fixed A and t 
we have 

pt-lim{2? A A -2? t , A } -0. (92) 

Since Z A X preserves the trace and is completely positive, the same holds true for T> A X in . This implies that 
II -^T A II = 1 ^ or an y ^ ( see ^ ne discussion after ([69])). Applying Lemma [TJ we deduce from (|92|) that 

2? A (T)^2? t , A (T) (93) 



for any T € Bi(Hp) with finite support on the lattice Z d . Hence 2\a can be extended to a bounded, trace- 
preserving and completely positive operator on Si ("Hp). From this, we straightforwardly deduce that, for a 
convergent sequence pp — > pp in Bi^Hp), 

lim V A x {p A ) = V t , x (p P ) . (94) 

AfZ d 

Since e -1 *^''! -> e -1 * 1 ^'" 1 strongly, (JM]) also holds if we replace X> A A by Z t A A = e" it[H p^ ] V A x and T> t ,\ by Z tjA 
which is given by (|7ip . Therefore, the limit in Proposition Q] exists and JE^a is trace-preserving and completely 
positive on Si (Hp). □ 

4.7 Crossing diagrams vanish in the van Hove limit 

In this subsection and the following ones, the correlation function f(x, t) do not need to satisfy the assumption 
(|32|) and we only require that 



||/0v)||i = / d*|/(x,*)| < oo for any a; G Z d . (95) 
Jo 

As announced in Step H of the plan of the proof, we show: 
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Proposition 4 Assume the integrability condition (|95p . If n is a crossing diagram, i.e., n 7^ 2Liadder> then for 
any fixed n, t, x, I, and ccq , j/o ^ x, y, 



lim / At \\V x , n (lL,t,x,l)xo,vo-:x,y\\ = ■ (96) 

A->0 J 

Proof. It is a simple adaptation of the arguments used in [35], Section 6.3. Let 7r = {(ti, a%), • • • , (i n , cr n )} be 
a crossing diagram. This means that one can find two pairs (i M , a^) € n and {i v , a u ) G 7r such that t M < t v (i.e., 
H < v) and < ct m . According to (j53"|). we need to show that 

I f m (t* m -Ol — > (97) 

Wa-»t) m=1 

as A — > 0, where we abbreviated f{x am -i lm ,t ffm — t tm ) by f m (tcj m —t Lm ) (recall that here x is fixed). One has 
J\{Z,T)< [] / di tm di ffm \ 2 \f m (U m -t lm )\ 

x / At h At^At L ,At aij A 4 |/ M (^ - t 4 „)| !/„(«„,, - 0| . 

The first (product of) integrals on the right-hand side is bounded by (sup m ||/ m || it)™ -2 . To deal with the 
last integral, we first bound the integral over t au by A 4 1 1 f v \ \ 1 / M {t„ — t t ) and then substitute v — \ 2 t a ^, 
w = A 2 (< ff)j —tt u )t an d t' = t a ^ — t Lfi . This gives the bound 

\\M\t [ Av [ Aw [ At'lf^t')] . (98) 

JO JO J\- 2 w 

A -2 

For any (v,w) £ [0,r] 2 such that < w < v, f^-a w dtf \f t i{t')\ converges to zero as A — > (because ||/ M ||i < 00). 
This integral is also bounded by H/^Hi, therefore (|98[) converges to zero by dominated convergence. □ 

4.8 Contribution of the ladder diagrams 

In this subsection, we determine the contribution of the ladder diagrams and accomplish Step III of the proof. 
We first introduce the following operators on B\(rlp): 

U x (t) = e - it[Hp '^ = c -^ 2 ^op+s,-] ( Uo (t) = e -n[S,-] _ (99) 

Define the family of operators 

/>oo 

A(z',l';x,l) = - / Atf(x'-x,t,l)U (-t)l(x',0,l')U (t)l(x,0,l) . (100) 
Jo 

The integral is convergent by the integrability condition (|95l) . Let 

<— 

W\,n(Z,X,l_)= Y[ Ux(-X- 2 T j )A{X2 j ,l2 j ;X2 j -l,l2j-l)Hx(X~ 2 T j ) (101) 
3=1) •••>« 

for any r = (n, . . . , r„) g E™ . We have 

Proposition 5 Assume the integrability condition (|95p . For any fixed x, t and t , 



lim 



/ dtV A ,n (bladder) t.i X, ,71 ^.5 

I) =0. (102) 



^2„(A- 2 r) Z„(r) 
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Proof. Let us set, for any 6 > 0, 
Ax 



Then 



/•A 

iS (x',l';x,l) = - / di/(x' - x,t,l)U x (-t)l{x',0,l')U x (t)l(x,0,l) . (103) 
Jo 

lim^ Ai 5(x',r;x,0 = .A(x', /'; x, Z) in norm and \\A x ,6(x' ,l';x,l)\\ < \\f(x' - x, OHiH^ 2 (104) 

A — >0 

for any x, x' <E Z d , and € {L, R}. This follows from the dominated convergence theorem, using (i) the 
integrability of |/(x' — x, - )|, (ii) the norm convergence lim A ->o (*) = Wo(i) (which follows directly from the 
boundedness of flhop), and (hi) the bounds ||I(x,0, Z)|| < ||W|| and ||ZA(i)|| = 1. Now, using X(x,t,l) = 
Ux{— t)I(x, 0, l)U\{t) and setting Sj = t 2 j — tsj-t, we rewrite the (infinite volume version of) (|62l) as 



(105) 



VA.nfeladde^X,!) =(-A 2 )" \\ f(x 2j ~ X 2 j-1, Sj, kj-l)Ux{-t 2j -l) 

j=l,...,n 

Ux(-Sj)l(x2j,0,l2j)Ux(Sj)l{x 2 j-X,0,l 2 j-x)U X (t2j-l) ■ 

We now perform the variable substitutions Tj — A t 2 j-i, Sj = t 2 j — f° r 3 = 1) * " , ^, to get 

/ dtV x , n (lL iaddc r,t,x,L)- J drW A ,„(r,x,{) = / dr { f[ U X {-^) (106) 

Z 2 „(A-2 r ) Z„(r) Z n (j) J =1 >-,« 

A\ tTj+1 -T j (X2j,hj;X2j-l,hj-l)U\&)- ]J W A(-^M(x2 J ,%;X 2 j_l,?2i-l)WA(^)} ■ 

j=l,...,n 

For any < t\ < t 2 < ■ ■ ■ < T n , the integrand inside the curly brackets converges in norm to zero because of 
(|104|) . This integrand is bounded by 2(max 7 -=i,... , n ||/(x2j — X2.7-1, Ollill^ll 2 )™- Hence an application of the 
dominated convergence theorem yields the result. □ 

4.9 Spectral averaging 

To end the proof of Theorem [TJ we use some standard techniques of "dynamical spectral averaging", originally 
used in Ref. [8] in the same context. 

Lemma 4 Let A and B be bounded operators on a Banach space y such that (e* B ) te H is a one-parameter group 
of isometries on y and the norm limit 

= lim r 1 f due- uB Ae uB (107) 



exists. Let £>(•), £*(•) be in C 1 (M,S(y)) (continuously differ •entiable B(y) -valued functions). Then, for any 

T>0, 

(1) lim/ dTiD(T 1 )e- ( - T1 ^ B Ae^ 1 ^ B E(n)= [ dn D{t x )A^E{ti) 
£ ^°Jo Jo 

(2) lim e -(r/s)B e (T/e)(B+sA) = Hm e (r/ e )(B+ e A) e -(r/e)B = q tA* 
e^O e^O 

Proof. To show the claim (1), we put E(t) = 1 (the general result follows by an obvious extension of the 
proof). Let us write D'(t) = £pD(r), then 



J dT 1 D(T 1 ) C - ( - Tl/e)B Ae^ T ^ E '> B = dn (d(0) + ' dr 2 D'(t 2 )\ e-^ Tl ' e ^ B Ae {T 

= f dr 1 D(0)e^ Tl / E)B J 4e( Tl / £ » B + f At 2 D'{t 2 ) f dn e^^Ae^ 

JO JO J T2 

£ Q (tD(0) + J dr 2 (r - T 2 )D'(T 2 )^j A* = J dn D(ti)A* . (108) 



(ri/e)B 
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To get the last line, we used (|107[) to estimate the integrals over n for all T2 < r, together with the dominated 
convergence theorem (since D'(-) is norm continuous and e tB is an isometry). 

The claim (2) can be proven from (l)by expanding the two first members as Dyson series in A (alternatively, 
see [36], Theorem 5.11, and the review by Derezihski and Fruboes in [20] for the same result under weaker 
conditions). □ 

Let us now apply Lemma [4] to the case at hand. We choose B = i[S, •], e = A 2 and y = Bi(Wp). Then 
(e tB ) t6 R is a group of isometries on Bi(Hp). Since S has a finite number of eigenvalues, the existence of the 
norm limit (|107l) is automatic. Moreover, for any A E B(Bi(Hp)) one has 

A^= V " AV " ( 109 ) 

where V u are the spectral projectors of [S, •], i.e., 

1 3 »{ T ) = SE a -E 3 „ u \s){s'\{s\T\s') , TgBi(Hp). (110) 

s,s'=l,-- ,N 

We first prove 

Proposition 6 Let us define 

!Cn(z,x,l)= I] U ^- T i)[ A ( x *hW,^j-iMj-i)^U^Ti) with U k (r) = e-^lor,,-] (m) 

3=1, ...,n 

where H^ op is given by L31\) . Then for the norm topology on B(Bi('Hp)), 

lim / dr W a ,„(t,z, f) = / dTK. n (T,x,l) 

A-s-0 J J 

Z„(r) Z n (r) 

Proof. An explicit calculation yields 



lim/ drWi \(t,x,1) = Urn / dn ZA(-Ti/A 2 )„4iW A (n/A 2 ) 

= lim / dr 1 ^(-r 1 )Wo(-ri/A 2 )^ 1 W (r 1 /A 2 )^(r 1 ) - f dn ^(-rO^W^n) 



(112) 



[ffhop.f = K op ,-]. (113) 
Choosing A = i[iJh op .-], the claim of Lemma 0] yields 

lim \\U q {-t/\ 2 )U x (t/\ 2 )-U^t)\\ = lim ||W a (t/A 2 )Wo(-t/A 2 )-^(t)|| = 0. (114) 

A — ^0 A — ^0 

We use the abbreviation Aj = A(x2j, hj', X2j-i, hj-i)- Since Uo(t) is an isometry, this gives 

hm o ||W A (-r/A 2 )^ A (r/A 2 ) -^(-t)Wo(-t/A 2 )^ (t/A 2 )^(t)|| - . (115) 

To prove the Proposition, we consider first the cases n = 1 and n — 2 and then conclude by induction. 
• n = 1: one has 



(116) 



where we used (|115[) in the second equality, relying also on dominated convergence and the uniform 
boundedness of the integrand, and claim (1) of Lemma|4]in the third equality (note that U^(-) is C 1 ). 
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lim / dTWa,A(l,£,i) = lim / dr 2 / ' d n U x (-t 2 / X 2 )A 2 U x (t 2 / X 2 )U x (~ Tl / X^A^n/ X 2 ) 



lim / dr 2 U x {-T 2 /X 2 )A 2 U x (r 2 /X 2 ) I ' dn ^(-n)^(n) 



A^O 



JO 



^2 



lim / dT 2 U^-T 2 )Uo(~T 2 /X 2 )A 2 U Q (T 2 /X 2 )U^T 2 ) / dri^(-Ti)i ^(n) 

A— >0 



f dr 2 w„ {-t 2 )A\u^ (t 2 ) / 2 dn W< (-Ti)^Wf, (n) 

Jo Jo 



(117) 



The second equality is the case n = 1, the third is (|115[) . and the fourth follows from the claim (1) of 
Lemma 2J 

• n > 2: follows by a similar induction step. 

□ 

End of the proof of Theorem [TJ Collecting Propositions [4] [5j and|H we have 

lT E / dt_(V x , n (K,t,x,l)) XOtyo ^ y = J dr(JC n (r,x,l)) XOtm (118) 

- 2 2 „(A- 2 r) Z„(t) 

for any fixed n, x, I, and x,y,xo,yo. Thanks to the second statement of Proposition [3] and to the dominated 
convergence theorem, we obtain in view of (|66p 



pt-limP A -2 T)A =^^ f dr/C n (r,ai,/) =e iT[H >»p'- 1 e^ ir[ ^op^ +T - Ah (119) 

where the last equality comes from a Dyson expansion in powers of A^ (the series is convergent in the pt-lim 
sense by dominated convergence), and 

A* = J2 [A{x,l;x',l')f . (120) 

x,x' ,l,V 

One concludes by invoking Lemma [T] (see also the discussion after this lemma) that 

V X -2 TA -> e iT[ ^ p.-] e -^Kop.-]+ T ^ strongly as A -> 0. (121) 
But c i *[' s ' ] and e _i *[ ffp ^] are isometries on Si (Hp). Therefore, in view of ([29]) and ([7I]l . 

PW(T) = lim e ^- 2 r[S,-] e -^- 2 r^H hop+ S,.] Vx _ 2T ^ = e -ir[^ op ,.]+r^ (pp) (122) 



in the trace-norm topology, where we have used again Lemma |4] (2). 

To establish the agreement with the generator & given by (|33|) . we check by inspection that, for any 

[A(x,l = L;y,l' = R)]\ PP ) = £ c(y - x, u)W u ® \x)(x\ PP W* ® (123) 

«eo-([-,s]) 

[.40M = £;y,r = L)]\pp) = W,, v Tp P + * x ,„ £ c(0,w)W w W* ® |x)<x| p P (124) 

uea([-,S]) 
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and 

[A(x,R;y,L)]\p P )=( y [A(x,L;y,R)]\p P )y , [A(x,R;y,R)]\p P ) = ([A(x,L;y,L)]\p P )Y . (125) 
Consequently, we have 

A\ PP ) = i[T, PP ] +Z(pp) - ^{2t*(l),p P } (126) 
with 21 as defined in (1361). □ 
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A Appendix: proof of Proposition [2] 

We show in this appendix that the main hypothesis (l3"2"j) of Theorem Q] concerning the decay of the bath 
correlation functions fi(x,t) is satisfied under the following conditions: 

(i) d > 2; 

(ii) the support of <7o,i(<?) belongs to the open ball B r = {q e T d ; \q\ < r} with < r < n; the form factor 
go,i(q) and the momentum occupation numbers Ci(<?) depend only on |g| on B r ; 

(iii) the bosons have a linear dispersion relation: Vi{q) — \q\ for q G B r ; 

(iv) the non-negative functions 

= IgoAvTCiiq) , ^,-(\Q\) = \9o,i(l)\ 2 ^ + Q(l)) (Al) 
belong to C 2 (]0,7r]) and the three functions of \q\ below are in L 1 ([0,7r]): 

i<?r n{d - 3 ^ } v>i >± (M) , \q\ d - 2 ^±(\i\) , i?i d -v-:±(ki) . (A2) 

For simplicity we omit the index i labelling the baths. Using the notation (|A1[) . the correlation function 
(|26|) reads 

f(x,t) = f+(x,t) + /_(*,*) = J (V+(|g|)e^e to M + ^-(kDe-^e-^W) . (A3) 

We first show that under conditions (i-iv), there exists a constant C<; > such that for any n £ N*, 

/ di sup |/(a;,t)|e-— < C d if rf > 3 , / dt sup |/(a;,t)|e"— < C^v^ if = 2. 

x£Z d ,\x\>t/2 Jl xeZ d .\x\>t/2 

(A4) 

Actually, by (|A3|) and (ii-iii), f(x,t) can be rewritten as the Fourier transform 

AM) = - ^ ^ /' da;|a;|^sign( w )(kl)^(|^|)e itw (A5) 

where sign(w) = ±1 for ±w > and J m (r) = (r/2)"T(m + i) _1 (2/\/7r) / Q du(l - u 2 ) m ~^ cos(ru) denotes 
the Bessel function of order m (here T is the Gamma function). A standard bound, see e.g. [37], yields 
C = sup^Q-tv 7 ?! J m (r)\} < oo. Hence 

sup |/0M)| < i ^ d-i / dw|w|V^ sign(w) (|a;|) . (A6) 

|aj|>*/2 227T2£ 2 J_ r 
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The last integral is convergent thanks to assumption (iv). Then (|A4=|) follows from 



f d 1 1 f d 1 

sup / dtt^^ e~^ < oo if d > 3 , sup -= / di i~ ^ e~ ^ < oo if d = 2. (A7) 

nGN* Jl nGN* V n J I 



We now show that there exists a constant C > such that 

C 

xez d ,\x\<t/2 " v ' " ~ t 



max ..,J/(M)I < 72 ■ ( A8 ) 



Let us set 

%,M) = Vz%) + f (A9) 

where V is the gradient with respect to g. Note that £, i) is the gradient of the phase S l (<7, a;, t) ~ i>(q) + q-x/t 
appearing in the oscillatory integrals (|A3|) . Assuming \x\ < t/2, one has \v(q,x,t)\ > 1/2 (since V^(q) = q/|g| 
has norm one), thus this phase has no stationary points. Noting that e ltS (i^ x ^) = (it) -1 (i>/|i>| 2 ) • ^ e ltS (i,x,t) anc j 
integrating twice by part yields 



(A10) 



where dk is the derivative with respect to the fcth component of q. Note that the boundary terms vanish thanks 
to condition (ii) . A simple calculation shows that the integrand in the right-hand side of (|A10[) is bounded by 

ci|<?rV±(M)+c 2 |gr 1 |VV'±(M)l+C3 max |a fc ^±(|ff|)| (All) 

for some constants ci,C2, and C3 > 0. The last function is integrable by assumption (iv) and this proves our 
claim (|A"8]) . 

Collecting the above results and recalling that t)| < f(0, 0), we conclude that 
1 f°° M 1 1 f°° |x| 1 f°° 

- dt sup |/(x,t)|e-— < -/(0,0) + - / dt sup |/(ac, *)|e — s~ + - / dt max \f(x,t)\ (A12) 

"-Jo n n Jl x,\x\>t/2 n Jl x,\x\<t/2 

converges to zero as n — > 00. This proves Proposition [2j □ 

Let us stress that conditions (i-iv) are not optimal for the hypothesis (l3"2"j) to hold. In particular, the rotation 
invariance of go and £ in (ii) and the linear dispersion (iii) have been chosen to simplify the proof and could be 
omitted at the expense of using the stationary phase method to evaluate the integral over the manifold v(q) = lo 
in (|A3|) . instead of using (|A5I) . In contrast, the first condition d > 2 is crucial: in dimension d = 1, (1321) is not 
fulfilled. To see this, let us assume that (ii-iii) hold and that the baths are initially at thermal equilibrium, i.e., 
= - I)" 1 . For d=l, (jMl) gives 

/r j iq(t+x) 1 „ig(t-x) 

_J^ 2 |e».-!| • < A13 » 

It is clear that sup^gz |/(af, t)| does not decay to zero at large times t. As a result, the integral in (|32|) diverges. 
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